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ABSTRACT 


This thesis is concerned with developing an optimal launching 
schedule for ABM's deployed for defense of a number of point targets, 
i.e. an ICBM complex. Let the attack occur in N stages with an 
offensive aeanees of saturation. A dynamic programming model is 
developed for formulating the problem and a linear programming model 
is used in its solution. Equations are developed for determining the 
optimal number of ABM's to launch on each stage of the attack so as to 
maximize the expected number of silos surviving. Multivalued point 
target defense is discussed and formulated but no specific solutions 
are offered. The ABM system (including radars, computers, etc.) is 
not considered subject to attack. Some discussion of this aspect of 
the problem is offered. Single ABM launchings per re-entry vehicle 
are assumed. A test is developed to determine for what parameters 


single launchings are preferred over multiple launchings, under the 


assumptions of the attack. 
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I. INTRODUCTION 


ed 


This thesis addresses~the problem of efficient allocation of anti- 
ballistic missiles (ABM's) during an enemy attack on United States’ 
land-based intercontinental ballistic missiles (ICBM's). The targets 
(ICBM silos) aeay pice as "point targets", i.e., the area of the 
target is small compared to the destructive capability of the attacking 
missiles. This paper hy concerned with optimizing BS Suchine strategy 
only, and reasonable assumptions are made concerning the strategy of 
the offense. - 

For a discussion of offensive and defensive strategy optimization” 


’ . 
< ~~ 


-during a single battle, the reader is referred to Perkins, F. M., 


ae 


> 


Optimum Weapon Deployment for Nuclear Attack [Ref 1] and McEwen, Ww. R., 
The Attack and Defense of Targets by Missiles [Ref 2].* Perkins assumes 


Feoae knowledge both by offensive and defensive forces, and a complex 
of point targets as the object of the battle, whereas McEwen discusses 
the problem of a single point target using both probabilistic and game 
theory models to obtain a solution. For a discussion of offensive 
strategy optimization only, refer to Piccariello, H. J., Missile 
Allocation [Ref 3] and McLaren, M. D., Walkup, P. W., A Missile 
Targeting Problem [Ref 4], and A Multiple Assignment Problem [Ref 5]. 
Piccariello takes into account the vulnerability of control centers 

as well as the silos. He develops solutions for the continuous and 
discrete cases. Both papers coauthored by McLaren and Walkup are 
closely related and should be read together. The first uses Monte 


Carlo techniques for estimating expected damage and minimum damage 


level for a given strategy. The second paper is concerned with optimal 
programming and tateociee ef missiles prior to an attacks 

This paper considers the attack to occur in stages, the number of 
stages being dependent upon the size of the missiie inventory of the 
offense. The question for the defense is, "What quantity of defensive 
forces (ABM's) should be expended at each stage of the attack in order 
to maximize the number (or value) of undestroyed silos after the attack?" 
Two general scenarios are considered. First, all silos (point targets) 
are considered to be equivalent in target value and second, all silos 
are not equivalent but rather belong to ordered classes. Both dynamic 
and linear programming methods were applied in an attempt to formulate 
and solve this problem. Although both of these methods were satisfactory 


in the formulation phase, they proved less efficient than a direct ana- 
lytical approach in obtaining a closed form oration 

From a practical viewpoint, the solution should necessarily be 
restricted to being inter-valued, since it would be meaningless to 
Toque the firing of 1/2, 1/3, or any other fractionseesan ABMwat aa 
attacking missile. This aspect of the problem is ignored but unlike 
Pennington's [Ref 6] solution to a similar problem, the integar con- 
dition 1s NoOLetnous@t FO be Criletcal “umes HERG. nacance of 
assumptions made concerning the geometry of the targets and the capa- 
bilities of the defense. Equations are developed, in terms of initial 
parameters, for the ee reales een of the ABM's. 

Interest was developed for this thesis by Pennington's paper [Ref 6] 
concerning the derense of a single point target against successive 


missile attacks. Unless a large number of ABM's were available for 


each point target, his solution often resulted in the firing of 


fractional ABM's. Although such a result is mathematically correct, it 
is not very useful. For example, firing 1/3 of an ABM at each of three 
attacking re-entry vehicles (RV's) results in a higher target surviva- 
bility than firing one ABM at only one of those three attacking RV's 
(see Appendix A for a proof of this statement). 

This analysis is restricted to single ABM launchings against any 
particular RV. If a sufficient number of ABM's is available to allow 
multiple launchings at a single RV, and to defend every undestroyed 
silo on every attack, then the equations developed herein are inapplic- 
able. An exact relationship defining the region of applicability for 
this solution is developed in section (III.B.8). The general method 
of this paper could be extended to a situation involving large numbers 
of ABM's if a solution to that problem were desired. However, the 
practicality demanded by scarce resources indicates that this would 
probably not be the case. It seems unlikely that an attack would occur 
at all if the defense were so well endowed with antiballistic missiles. 
This implies a rational behavior on the part of the offense, which 


might or might not be justifiable ina real situation. 


IT. PROBLEM DESCRIPTION 


A. ASSUMPTIONS 

If and when an attack should ever occur against the United States' 
intercontinental ballistic missiles, and should an antiballistic 
missile apes be installed to defend those ICBM's, there arises the 
obvious question, "How does the defense allocate ABM's so as to pro- 
tect as many silos as possible from destruction?" In order to answer 
such a question it is necessary to make some assumptions about the 
nature of the ne elen the reliability of the equipment, and the 
APS (eos EE the situation. 

A "farm" will be defined to be an area of land throughout which 
are scattered a number of silos containing ICBM'’s. The separation 
of the silos is sufficient to disallow aanlieiterie destruction by a 
single oe vehicle. If the re-entry vehicle is sufficiently 
close to one silo to destroy it, then it is too far from any other 
Soe ao cause its destruction. A re-entry vehicle is a single bomb 
and rhemeaee that it might have previously been released from a 
warhead containing many bombs is not significant. Technology 
available to distinguish decoys and similar systems from actual 
re-entry vehicles, as well as the space and weight demanded by such 
systems preclude their consideration. Every radar target approaching 
the farm on a proper trajectory shall be treated as a re-entry vehicle. 

The attack is to occur in stages, with the offense attacking every 
silo in the farm on every stage. The offense will not be allowed a 


"shoot-look-shoot" capability, i.e., the ability to ascertain battle 


damage between stages and adjust his strategy accordingly. Since each 
silo has the same probability of being undestroyed after the first 
stage, the offense is necessarily required to attack all the silos 
again on the second stage, or none. Unless the silos are of different 
strategic value (to be discussed in chapter IV) they are equivalent 

in the attacker's eyes on each stage of the attack. One further 
assumption ig Kieties this reasoning more completely. Each ABM in 

the defensive system has the capability of defending any silo within 
the farm. Thus it is not possible for the offense to concentrate 
their attack so as to exhaust defenses in a portion of the farm. If 
there is any ABM left for defense, it can be launched to intercept an 
RV aimed at any — within the farm. Thus the offense, if it attacks 
the farm at all, is assumed to attack each silo at each stage until its 
inventory is exhausted. 

The silos, and the ICBM's they house, are the only targets attacked 
by the offense. All ICBM control centers, the ABM complex itself, and 
all associated radars are eliminated from consideration as targets. 
Chapter IV, part B, discusses this subject more completely. 

It is assumed that the defense can ascertain battle damage between 
stages and thus defend only those silos which remain undestroyed. Once 
a silo is destroyed, it is left undefended on future stages of the 
attack, thus all RV's aimed at destroyed silos go unchallenged by the 
defense. If a RV successfully penetrates the defense but misses its 
preassigned target, it misses all targets. No lucky hits are allowed. 
In a similar fashion, it is assumed that one ABM can destroy only one 


RV. The spacing between offensive stages and between RV's within a 


stage is sufficient to disallow multiple successes by any ABM. The 


spacing between stages is not sufficient, however, to allow launchings 
of ICBM's during the attack. Only those ICBM's which survive the 


entire attack can be used by the defense ina retaliatory attack. 


B. DEFINITIONS 
pete. terms as they are used in this paper have been defined in 
part A of this chapter. However, for ease of reference and complete- 
ness, all terms used which require an exact definition are listed 
below. 
Point target: a target which has an overall surface area which 
is small when compared to the destructive capa- 


bility of a single re-entry vehicle. 


ABM: antiballistic missile. 

ICBM: : intercontinental ballistic missile. 

RV: re-entry vehicle, a single fission or fusion bomb. 
Stage: a single wave of RV's which all arrive at their 


targets at approximately the same time. There is 
one RV per target per wave. 

100% defense: that defensive strategy such that all undestroyed 
Silos are defended on a given PP 

Attack: a sequence of N independent stages, with stage N 
occurring first in time and stage one occurring 
at the end of the battle. 

Teme, ° the expected number of undestroyed silos immediately 
prior to stage (N-i) of the attack. 

rr the number of ABM's remaining in inventory immedi- 

ately prior to stage (N-i) of the attack. 

the probability that an ABM launched on stage (N-i) 


successfully intercepts and destroys an RV. 
10 


qc the probability that an undestroyed or unchallenged 
RV destroys its assigned target. 

The fact that stage N of the attack is defined to occur first in 
time and stage one to occur at the end of the attack might seem strange 
at this point. The reason for this backward numbering comes from the 
dynamic programming formulation of the problem, and rather than use a 
different notation in other sections, it is used consistent ly gor val | 


formulations, and for the solution. 


C. A DYNAMIC PROGRAMMING FORMULATION 
1. Notation 

At every stage of the attack, there is an expected number of 
silos destroyed. This destruction can occur for a number of reasons. 
First, every RV may not be challenged due to a scarcity of ABM's. 
Second, except as a limiting case, ABM reliability for successful 
interception Nomicas than one. Of the RV's that go unchallenged and 
those that successfully penetrate the defenses, some will be success- 
ful in destroying their targets. Let Yea (stage return) be defined 


as the expected number of silos destroyed on stage (N-i). A diagram 


of the entire attack is shown below in Fig. 2.1. 

















d 
N eH Be 
X X Xy ! a 
N . N- = -i . a (it) s 
T T aE oe) ae 
N N-1 = N-i ss N- (itl) a 
y | | 
Y 
N YNed Y 


ting. 


Figure 2.1 
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Each stage of the attack is identified by the index (N-i), i = 0, 
1, ..., N-l. At the beginning of each stage, there are a number of 
ABM's 'st PIM in invencory, Ryu and an expected number of silos still 


undestroyed, T < For this stage, a quantity of ABM's, d 


N ~~ is 


N-i 


launched, leaving XN (41) ABM's for future defense. Yy-; Sites are 


destroyed (expected number) leaving T as the expected number 


N- (it+1) 
of undestroyed silos for the beginning of the next stage. The d's 
are the decision variables which are to be selected in an "optimal" 
manner. 
2. The Mathematical Model 
The problem, as stated in these terms, is then to select the 
vector .d = (d,> d,, sae di) such that the sum of the expected number 
of silos lost on each stage is minimized, i.e. 
N 
mitt Deu by selecting €2.1) 
qd 
J 
d d oageees di in an optimal way. 
a. The Stage Transformation Equations 
Since dyed ABM'’s are launched on stage (N-i), the number 


remaining at stage (N-i)-l is 


#5 1a) ee ae (2225 


and also since ty i silos are expected to be destroyed on stage 


(N-i), the expected number remaining is 


= = o Ze 
ty i+1) eke Sea: - (2,3 ) 
However, Y, . is a linear function of T, ., and d, . as will be shown 
N-i N-i N-i 
below. 
(1) An Expression for Yea’ Of the dad ABM's launched 


in defense of stage (N-i), are expected to be successful in 


PN-i 


12 


destroying their assigned RV's. But Ty RV's are approaching the farm, 


of which Tea are directed at undestroyed silos. Since the defenders 


are concerned only with defense of undestroyed silos, the ABM's 


launched are directed at those Ted RV's aimed at the undestroyed 


silos, and the other RV's are left unchallenged. Thus, of the Ty-i 


RV's that are challenged T d are expected to penetrate 


N-i ~ Py-i °N-i 


the defenses. Of these, the probability that each destroys its 


~ 


assigned silo is q and thus, the expected number of ICBM's destroyed 


on stage (N-i) is 


Yyeg ~ DO ys 7 Pyeg Sn-a?- i as 
Therefore 
ides her ~ ME “See Se? 
or 
Ty- itt) * Tyg 7) + oPyig Sper ee) 
3. Analysis 
Let £,(X,, T,) = min (,) (2.6) 
d 
i 
where Xi» Tt) = de > 0 
and £, (X, T,) = min (Y, + £1 (X,; T,)) MT 
d 
vi 
where Xo » qT, 2 d, 2 0 
In general, let 
f. (X., T.) = min (¥. + £. me ps A. 2.8 
j Bye By win + Oy ys TDD (2.8) 
J 
where X., 1, 26cm >: 
J (nh | oe 


13 


This expression, (2.8), along with the stage transformation 
equations; (2.2) and (2.5) is a dynamic programming formulation of 
the problem. Unfortunately, the problem cannot be solved in closed 
form as presented here. To realize this, consider eq. (2.6) 


f) (X,; T,) = min (Y1) 
a 


tt 


fi (X5 Tt) min (qT, - qp,d,). 


a 


ad 


Since the coefficient of dy is always non-positive, and negative for 
all practical cases, the quantity in parenthesis can be minimized by 


making d, as large as possible. However, d, cannot be larger than 


1 1 


X, since no more ABM's can be launched than there are in inventory, 


and d, cannot be larger than T, since it is unnecessary to defend 


1 1 


any Silos which were previously destroyed. Now the matter becomes 


a question of whether X, > Ty or tT 2 Xi. 


Consider Xx) > T,- In this case there would be more ABM's 


available for defense than were needed (recall only single ABM 
launchings are allowed against any given RV). This implies that 
the excess ABM's should have been used at an earlier stage as 
Opposed to not using them at all. Any strategy, to be optimal, 
will use all available ABM's for defense. Thus X, must be equal 


i 


Eo or less than T, and 


where d¥ is that value of d, which satisfies £, (> T,)- 


Continuing with eq. (2.7) 


f, (X, » T,) = min (Y, + f, (X1> T,)) 
d 
2 
= i i * 
f, (xX, , T,) a (Yo + qT, qp,d*) (2.9) 
Z 
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where Xo» T, 2 d., 


Equation (2.5) with i = N-2 is 


> 6. 


and eq. (2.4) with i = N-2 is 
Yo = q(T, - Podo)- 
Then eq. (2.9) becomes 


2 
— j ~ - = = * 
f, (X,, T,) = min (a2 q) T, - (apy - 4 podd, - 4p, ). 
2 


since 1 > q > 0, thetceefficient of d, is always non-positive, and d., 
must assume its~maximum value in order to minimize the expression in 


parenthesis. This is where the problem arises if a closed form 


solution is desired. It is not known whether X, > T, or T, > X 


2 iz 2 2 


and since the smaller of these two quantities is the upper bound 
FOr d,s the eae value of d, cannot be determined immediately. 
A tabular solution could be obtained once initial parameters were 
specified but a more direct solution is sought. 

The following chapter gives two linear programming models of 
this problem. The second model is solved in closed form for all 


decision variables. 


15 


TTI. LYNEAR PROGRAMMING FORMULATIONS 


The problem addressed in this chapter is exactly the same as 
introduced earlier. The same assumptions regarding the structure 
of the aeeaen and the defense will be used. All previously introduced 
notation is also unchanged, and additional notation is defined as 


required. 


A. MODEL I 
This model is presented for continuity and will not be solved in 
closed form. Therefore only necessary equations are shown and the 
details of their development are contained in Appendix B. However, 
the reader should be cautioned that there are many aspects of model 
I which are exactly the same as model II and an understanding of 
model T is necessary before going on to model II. 
Let 2 represent the expected total number of silos lost after all 
N stages have occurred. Then 
N 
A dae 


j=l 


tn 


ne (see (2.1)). (1) 


However, each is a linear function of He and i as shown by eq. 
(2.4). It can further be shown (see Lemma 1) that each tT, Tse) 
linear function of d's, i = jtl, ..., N; thus, the expression for 
z above can be written | 


Z = D, + Did SF Woke Duda (see Appendix B). (332)) 


Then the problem can be stated as follows 


minz =D +D,d, +... + Dy day (39) 
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subject to 


N-1l 
OY gest TON 

i=0 
(2) dag < Tea i) , N-l < (344) 
@) dua Zul) i=) OF euue sae Ie. 


The objective function, z, is linear. The coefficients of the 


variables (derived in Appendix B) are as follows 


a 


-| ’ 
pat, ( £ (iy!) 


He 


Ds = qp, ¢-1) 


e _ 7 J : 
Dy-i 7 Wy. (4 = (1-a)) = 1) (3.5) 


, N-2 ; 
Dy = Py (2 cea) ~ 1). 


1. Interpretation of Constraints 
Genreraint (1) requires the defense to launch all the ABM's 
in inventory and available. Constraint (2) is discussed in detail 
below and constraint (3) is a set of N non-negativity restrictions. 
a. A Modification of Constraint (2) 

Recall that the objective function is based upon the 
assumption that only one ABM is launched at any given RV. Constraint 
(2) aeeeres cat this assumption is not violated, but it must be 
modified to be useful since each oe is a linear function of some of 


d Constraint (2) in the 


the decision variables, namely rie veey Aye 


form used above is understandable in relation to the physical 


17 


Situation, but the form given below would have to be used if model I 
were to be solved using given parameters. 
N 


ay ee p.d,(1-q 
jen-iti J J 


pe aes) (3.6) 


ace I= 0,1, 7.4.5. 
This form of constraint (2) is stated as Lemma 2 in the next part of 
this chapter and a proof is given. 
If one wished to do so, this model could be solved by means of a; 
simplex algorithm. Naturally, all parameters must be specified. 
However, a general closed form solution for an equivalent model is 


presented in part B which offers a faster, more efficient solution. 


B. MODEL II 

1. Preliminaries 

The two previous models used to formulate this allocation 

problem both utilized the concept of minimizing the expected losses 
or sum of expected losses for the N stages of the attack. Consider 
now: an eqnivalent: model which is developed in order to maximize the 
expected.number of undestroyed silos at the end of the attack. Since 
the. defenses’ ICBM's are designed for retaliatory, or second strike 
capabilities, is important to have as many as possible remaining 
to carry ,ut thatemission. 

Recall that TS is equal to the expected number of undestroyed 
silos after the farm has been subjected to an attack of.N stages. 
The number of undestroyed silos prior to stage N is just Ty since 
none have yet been lost. It has been assumed that the offense attacks 


every silo on every stage, so the number of RV's aimed at undestroyed 


Silos on stage N (first in time) is also equal to Ty 
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2. Development of the Objective Function 
Using eq. (2.5) with i=l produces 
ee ee 
and for i=0 yields 
Tye] = Ty (l-q) + qPxy di. (3.8) 
Substituting eq. (3.8) into eq. (3.7) for Tye produces 
Ty-2 = Ty @eae + q(l-q) Prt = WP yFy td: (3.9) 
Comparing eq. (3.9) with eq. (3.8) leads to the general formula of 
Lemma l. 
Lemma 1 If there are Ty silos to be defended at the begin- 
ning of an attack, then the expected number of silos remaining at 
the beginning of stage (N-i), i=l, ..., N-1; is 
N 
Ted = Te(1=q)" +q 2 (19g) OY a, (Seo) 
jaN-itl 
Proof Lemma 1 will be proved by induction. Note that eq. 
(3.8) is exactly the same as the equation produced by eq. (3.10) for 
i=l. Therefore Lemma 1 is valid for i=l. Assume Lemma 1 is valid 


for i=k. From chapter II, eq. (2.5) gives 
Tye ett) ~ Tyan Sb @ + oPy iS: Con) 


Suibstitutimg eq. (3.10), using i=k, into eq. (3.11) for Tek yields 
k x j- (N-ke+1) 
Ty- (et) * Lbyft-a) + =< pd, (l-q) + apy _pdy_y: 


Combining terms produces 


N 
7 k+l j- (N-k) 
TNH (ek 1) Ty (1-4) ts ijn? P54 
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which is eq. (3.10) for i=ktl, and which completes the proof for Lemma 
il 
Therefore, the expression for TA which is the new objective 
function, is 
N : 
=: N gi-t 
max T = T_ (li-q) +q 2 (i-q) pa.. 3.125 
O N ‘ IP 3 
j=l 
3. Constraints 
The constraints applicable here in model II are exactly the 
same as those given in model I, and will be restated following Lemma 
2. Constraint (2) must be modified however, as mentioned in part A 
of this chapter. 
Lemma 2 For each i; i=0, 1, ..., N-1; the following inequality 


is a restatement of constraint (2) 


| N : , , 
dy 7a (a)? ma T,(1-4)°. (3.13) 
*  4=N-i+1 bee 
Proof Constraint (2) requires that, in general, 


Oe Sel ig 89 (3.14) 


However, from Lemma 1 


N 


Teg ss Glas b +q 2 (r-faOEe ow V, 65s No 1. 
“ jsN-i+1 JJ , 


Substituting this expression into eq. (3.14) yields 


N 


dy, < Ty(1-a)" tq Coo ec 
j=N-itl a 
Rearranging terms 
Bs {- (N21) i 
abo} -q 22 (1-q)3 p.d, < TyC-@) 
j=N-it1 : 


which completes the proof of Lemma 2. 
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4. The Complete Form for Model II 


In summary, model II is 


N 


max T, = Ty(1-a)" + qd acct 
sal saad 
subject to 
Nl 
ae pes 
1=0 
2 jai) i 
OG) GL re ¢ AEG ap pd, < T,(1-a) (3.16) 
j=N-itl 
i=0, ...aehel 
a > = - 
C3) deg > 0 TO , sees N= 


It is this model from which the closed form solution is obtained. 
De The Solution 
a. A Physical Consideration of the Objective Function 
The objective function can be written more concisely as 
N N 
T =T.(1- + 2% C.p.d, oe, li 
: reaalga ath F584? 544 Gri) 
J 
= ie N , 
where 2. = q(l-q) . Note that Ty (1-4) is a constant and has no 


effect on the optimization process once T, and q are specified. 


N 
Consider the case where 


C > C5Po ae CPx" (3.18) 


lat 
This inequality will be valid whenever 


* 


eT ; | 
(1-q)1 P > (1-q)/ Pia HOE J Lin ess 5 Nm 1 


ieeeriue, Of Pi/P aay > (l-q). 
The quantity Ps/P aay is a fraction between zero and one since it has 


previously been assumed that ABM reliability decreases in time. Also, 
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practical values for q should be in the neighborhood of 0.6 to 1.0. 
Thus, if p deteriorated by no more than 60% per stage from its 
previous value, eq. (3.18) is valid. It is hoped that ABM relia- 
bility would not be reduced so drastically as to invalidate eq. 
(3.18). 
b. A Characterization of the Optimal Solution 

Theorem 1 If xy ABM's are available a defense of an 
N-stage attack, and if as > C.P. for 1 < Wie t2ar2 > e.., N; then 

d¥ = 0 fore l—Kktl. 2.5, N 

d* > 0 


and d* = T. for i=k-l9ae.., lL for some N< k < l 
where % q* = xy and d* is the optimal value of d.. 
i Th 


Bio0 faeG is the largest coefficient in the objective 


eal 
function and TS can be increased most by making dy as large as possible 
since the coefficients of the decision variables in the first constraint 
are all equal to one. Likewise, CP5 is the next largest coefficient 

and resources should be allocated to d, whenever dy is at a constrained 


maximum. This reasoning can be continued up through C which com- 


NPN 
pletes the proof of Theorem l. 

In concept then, the solution is straightforward so long as the 
coefficients are in descending order. However, obtaining equations 
for the decision variables in order to allow direct computation is a 
‘more involved process. 
ec. Additional Notation 


let W represent the expected number of ABM's required 


to defend all undestroyed silos from stage n to stage one given there 


Z2 


are qT undestroyed silos at the beginning of stage n. To can have 
any value from zero to Tye 
Let Gs represent the expected number of ABM's required to defend 
all undestroyed silos from stage n to stage one given no ABM's were 
launched prior to stage n. 
d. Analysis 


From the definitions above, it is obvious that 


W, =T (3,.089)) 


Suppose xy — Wi Then the defense can achieve 100% defense on stage 
one. Next, suppose xy > Wie In this case, those ABM's not required 
on stage one are allocated to stage two. However, as ABM's are 


allocated to stage two, the upper bound on d, increases (see con- 


1 


straint (2)) due to the increase in T,- 


Let oa represent the expected number of undestroyed silos 
immediately prior to stage i given defense started on stage j, 


j >i. Then for 100% defense on both stages one and two 


Wo = qT, + T15 (3.20) 


But Ti5 is related to T, in the manner of eq. (2.5), specifically 


Tio = T, (1-q) + dPod,- (221) 
Substituting eq. (3.21) into ,eq. (3.20) yields 


Woes Tfl + (i-q)) + ap,¢,. 


However, from the definition of Wo (100% defense), d, in this instance 


equals T Thus 


a 
Wy = T, [1+ (1-q + ap,)]. (Ean) 


In similar fashion, the expression for W, can be written 


3 
W, = T, [1 + (1-4 + apg) + (1-4 + qp3) (1-4 + dpo)]. (G.23) 
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Lemma 3 For any stage j, N< j < 1, 


es 
te U1 ‘Sa emia Glee MI » (3.24) 
— i=1 k=l laun 


Proof For j=l, eq. (3.24) reduces to 


which is eq. (3.19). Thus Lemma 3 is valid for j=l. 


It is assumed that Lemma 3 is valid for j=n. Consider stage 


W = T +a T 


+ 
ot illest syepic i 


eg @n25) 


But fi is a particular value for Ti: Therefore, from Lemma 3 


orl 


the number of ABM's required to defend 100% from stage n to stage 


one is 
n-l i 
ss a ee G a ia ie Lal Pye) 
| i=1 k=l 
Also 
Sec al Fj tal Jacl qP rd Sasi © co) 
36) 
n-l i 
an ~ (Ty A-o ‘4 Pas itneL) G 7 q Ces Pate) eee 
i=l k=l 
Again from the definition of Wiel 
ri ;; Gosre 
and 
n-l i 
ee een G i" = —— a WP tet) Cnn 


Returning to eq. (3.25) 
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n-l i 


i 
| 


n+l Dey GQ eet aaa a Ta giommeer 


i=l k=l 
i 


| n 
ntl a G i 2 Wy Seal =" onan ae). 
i=l k=l 


= 
| 


) 
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which is equivalent to eq. (3.24) with j = ntl. This completes the 


proof of Lemma 3. 


Using the definition of a; and Lemma 3 produces 


ee G) Ty(l-a)y 4 
j 


since: T, = Ty(l-g) if no ABM's are launched on stages ntl, 


Using the first constraint and Lemma 3 with j=n-l yields 
n-2 i 
xy = dis ly G + 2 Tt (l-q + aP,-),): 
i=l k=l 


IPeeAe be the Coekricient of TO above. Then 


-1 
Xy =" gtr A. 
n n-lL 

Also 

Tiel = T (tea) + qpd and 

Ne-n 

T= TyG-@ 

Since no defense occurred on stages mtl, eee mS 


N-nt+l 
Xy dq, + (4-4) + apd) A- 


Solving for d* gives 


x. - AT Acc pier 
a=. = N N 
n [bap AqP 
Continuing 
N-ntl 
* = = 2 ' 
Sai cat ries ) , Ae 
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(i350) 


(3732) 


(3433) 


Similarly 


a = dx e 
dfugrs a*t, (l-aaewi ape | 2 (33343 
In general 
* fairs - va. a 


This form of the solution is most practical for actual computations 
since the equations for each ai in terms of initial parameters 
become very large: and unwieldly. 

e. Solution Algorithm 

In summary, the solution ale eerie is 

(1) Compute Qu? Que -. Q5 in that order using eq. (3.30) until 
that pair of Q's which bracket Xy is found. This yields the value” 
Guin 
(2) Compute d* ei eq. Ge3sF 


(3) Compute ds using eq. (3.33). 


lt 
(4) Pompmee ae j=2,. oe n-1; us ing eq. (3555). 
In a practical sense, it is only necessary to compute n and then 
d* and defend all undestroyed silos on each succeeding stage. Based 
on expected values, this allocation will produce the maximum surviva- 
bility of the silos. 
f. The Optimal Value of the Objective Function 
The objective function, eq. (3.12), was 
max T, z Ty (1-4) +q x (1-q)?""p,4,. 
j=l 

Because the optimal solution indicates that no ABM's are launched 
prior to stage n, all decision variables with subscripts greater 


than n are zero. Thus, the optimal value for the objective function 


can be expressed as 
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n ; 
| i 8 Gea” eq 2 (1-q))" ‘pa. (3.36) 
O N j=1 Jie 


However, from eq. (3.33) 


ae = TC-a) 4 gp ae (3.37) 


ne 


and from eq. (3.34) 


* = as = 
% d* (1 i 4P_1? 
or 
N-ntl1l 
ae = = * os 
d* (1 qd) + qp d*) (1-q + qp_j)- 
Similarly 


N-ntl 
* = = = = * 
d*_, (l-q tap, _») (l-q + qP,-v (a q) at ap,d* ) 


Lemma 4 For every j; j=l, ..., n-l; 
i 
N-nt+1 *) 
(leq + ap,_,) (ty(-a) Np a). (3.38) 
Proof Lemma 4 will be proved by induction. For j=l, 
@Ga5 (1.39) reduces to 
N-n+l 
23 = = * 
which is the same as eq. (3.37). Thus Lemma 4 is valid for j=l. 
Assume that Lemma 4 is valid for j=m. Consider d* From 


n- (mrtl)° 
eq. (3.35) 


yy = * = 
SHEE See bi i i 
But from Lemma 4 


m- 1 


A N-ntl *) 
¥ = - = eee). 
daa ee q+ qp._,) GyG-ao + qp.d 
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Thus 


m- 1 
N-nt1 
fe = = a om 7 
a¥ Get) 7 1-4 + QP yg) T I-d + 4P gg) (ty G-ay 7"! + ap ax). 
This can be written 
- N-nt1 
~*~ = = = * 
oa Gel) eg Cah te, (1,1 q) gi ap,a*) 


which is the same form as Lemma 4 where eae This completes the 
proof of Lemma 4. 

Returning to eq. (3.36) and substituting eq. (3.38) for each ae 
j=l, ..., n-l; yields 


N n-l 
* = = = * 
Ts Ty (1 q) + q(l-q) pas + BC (3.39) 


where B and C are defined as follows 


ie. n-(j+1) 


Eye 2 
_ q) P; Cr CE ee ape) 


ws) 
ll 


N-ntl 
2 te 
Ty(t q) + qp_d*. 


oe 
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Using eq. (3.32) and substituting for ds yields, after @actoring 


N-ntl 
-AT,.(1-q) : 
N-n+ n-1 - ) Ny N C neal 
k= — . ee = 2 
Te = Ty(1-q)™{ (1-q)"* + B) + ie ap, (Bt(1-q)""* 
(3.40) 
Equation (3.40) is the value of Lx as expressed in terms of initial 
parameters and n. This is a cumbersome form for calculations, but it 
is developed here for use ina later section. See section 7 fora 


more direct means of calculating T. 
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6. Discussion 
a. The Use of Expected Values 
With the aid of computers, it would be possible and indeed 
worthwhile, to recompute an entirely new solution following each stage. 
This would result ina new d* after each stage, and avoid the inherent 
discrepancies resulting from an analysis based on expected values. 

It is unlikely that the actual number of undestroyed silos present 
at any stage of a real attack would be the same as the expected number. 
It must be remembered that the parameters q and N are estimates based 
on the latest intelligence information available and are subject to 
error. Even the P; values, for which some data can be gathered, are 
subject to error. ABM's have never been launched or radars tested 
under conditions of severe atmospheric disturbance as would be caused 
by numerous nuclear blasts. Therefore, if an updated solution based 
on actual values could be obtained between stages, better results 
could be expected than by rigidly applying the complete initial 
solution. 

Consider the situation in which an updated solution cannot be 
obtained between stages. Based on the foregoing analysis, the 
solution that should be used would be to defend on stage n with d* 
ABM's. On each succeeding stage, defend all undestroyed silos, 
whether or not the number was greater or less than expected for that 
stage. If n were large, it is reasonable to assume that the number of 
stages requiring more ABM's than expected would have a cancelling effect 
on those stages requiring fewer ABM's than expected, and the shortage 
or excess of ABM's at the end of the attack would likely be small com- 


pared to X However, if n were small the excess or shortage could be 


N° 
a sizeable proportion of Xv 
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b. An Intuitively Appealing, but Erroneous Solution 
When this problem was first considered, one apparent 


solution was to select a subset of T, and defend that subset completely 


N 
throughout the attack. Since it was assumed that a sufficient number of 
ABM's would not normally be available to defend'all the silos properly, 
it was reasonable to assume that the number of ABM's available would be 


sufficient to defend a smaller portion of T In general, of course, 


N° 
this is not true. However, if q equals one, this is precisely the 

solution one does obtain from model II. So for this limiting value of 
q, ‘the intuitive solution is also correct. It would be foolish not to 


defend on stage N since every silo undefended in this case is destroyed. 


For q=l 


Therefore 
Qn 2 Ay = Sy-1 
for all eenisean it values of Xe Thus, defense starts on stage N in 
all cases when q equals one. 
c. The Peiiatt of N 

The estimate of the parameter N is critical in this 
solution. It is noc ne eant for the defense to determine N, and for a 
proper determination, the defense should be aware of the risks involved 
in overestimating and indaracuinaelte the parameter. If the estimated 
value for N tf: out to be greater than the actual value, ABM's will 
be left over, thus vaccine resources and losing more silos than 
necessary. On the other hand, if the estimated value for N is smaller 
than the actual value, the offense can attack unchallenged resulting in 


excessive losses of ICBM's. 
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Consider the case where q equals one. Underestimating N would be 
disasterous since one unchallenged stage results in complete loss of 
all remaining silos. Overestimating, on the other hand, would result 
in unused AMB's, but there would also be an expected number of undes- 
troyed silos remaining. Conversely, for q close to zero, the opposite 
is true, everything else remaining unchanged. Other parameters also 
affect this balance. Without specific parameters to consider, it is 
difficult to make any quantitative statements about the erroneous 
estimation of N. Decision Theory and Game Theory both offer techniques 
for optimizing the choice of N. A great deal of additional study could 
be done on this facet of the problem. 

7. A Numerical Example 

Consider a four stage attack and let the initial parameters 
of the problem be as follows. 


N 


a 


340 ABM's 


150 silos (ICBM's) 


N 

q = 0.8 
p, = 0.9 
Pp, = Owed 
Py = Oxo 
Pp, = O:..6%. 


a. Determining Q Values 


From eq. (3.30) 


4-4 


or Q), = Whe 


ot 


Substituting CF W, using eq. %@3.24) 
8 i 
= l Ee 
Q = T, Ci + =m (1 - 0.8 + Hyped) 
i=l k=1 


Using given parameter values 


Q, = 150 (1 + 0.92 + 0.74 + 0.53) 


or Q, = 477.9. 
Then 

Qa), ~ 478 > 340. 
Similarly 


. ee 
Q, = 79.7. 


=~ 80 < 340. 


op) 
O 
ye) 
WW 
Q 


b. Determining n and dé 


Since Q), = x 7 Q3> n=4, and from eq. (3.32) 


Be a 

340 - 150(0.2)( 1 +3 EH @We4cteOs8p )) 

< = al 4-k 
i=l k=l 
d* = 
n Z i 
1+ (0.8) (0.9) (1 + > T (0.2 + 0.8p,_,)) 
i=l k=l 


d*¥ = 99.1 = 99. 


n 
c. Completing the d's 


From eq. (3.33) 


dz = T, (0.2) + (0.38) (0. 9G) 
ds = Ol. — 201. 


From eq. (3.35) 


* = d*¥ = 
dy = d¥ (1 - 0.8 + 0.8p,) 
dy = 101(0.2 + 0.6) 
dF = 80.8 ~ 81. 


a2 


And also from eq. (3.35) 
di = 58.4 > 59. 
4 
Note that %& dy = 340. 
i=l 
d. Expected Number of Silos Remaining 
Using eq. (3.12), the objective function, and the results 
above, the maximum expected value for Ty can be determined. 
4 


150(0.2)° ~ Ons 3 (0.2)37+ 5 ax 
jel | 


max T 
fe) 


max TS SOMA Sii8\. 3 2ucte8 a8. 2 sheet 0.57 


Il 


max Ty 39.98 = 40. 

In summary, the solution above instructs the defense to launch 99 
ABM's on stage four, 101 on stage three, 81 on stage two, and 59 on the 
final stage. This allocation will result in saving approximately 40 
silos from destruction. 

As a sidelight, notice that the expected number of silos saved 
from destruction may be computed more directly without the use of the 
original objective function. Because of the nature of the solution, 
the following method is valid. Since d¥ = 59 ABM's, there are 
expected to be 59 undestroyed silos at the beginning of the last stage. 
Also, Pi = 0.6, so 35.4 or approximately 36 RV's are destroyed on 
stage one. Thus 23 RV's are expected to successfully pass through 
the defense, but of these only 80% are expected to be successful in 
destroying their targets. Therefore 35.4 plus 4.6 silos or 40 silos, 


are not destroyed on stage one. This is the same result one would 


obtain if eq. (3.35) were extended to include j=n, with dx = Tx, d* 


oie: 


has no physical meaning, but the pattern of the solution allows for 
this simple extension. Thus 


* = d*¥ = 
Ts d¥(1 eat qp,) 


which algebraically carries out the same operations which were performed 
verbally above. 
8. Maximum Number of ABM's Allowable 

Recall that the models presented in this paper are all con- 
cerned with single ABM launchings. Therefore, there exists an upper 
bound on the number of ABM's which can be used in any anticipated 
attack. Should the defense have a greater number of ABM's than this 
upper bound, multiple launchings would very likely be preferred over 
single launchings and the solution previously presented cannot be 
rire ta by applied. It remains to determine this upper bound. 

a. Development 

The maximum number of ABM's the defense can launch, using 

single launchings only, is that number for which every stage is 
defended completely, i.e. every RV aimed at an undestroyed silo is 
challenged by one ABM on every stage of the attack. Thus 


d* = T 


N 'N 

ane Coogee ~ Wy 

since 100% defense occurs on each stage following stage N because of 
the nature of the optimal solution. From Lemma 3 

bt: 


N- 
=T G+z 1 (a+ apy ,a))): 
Nee N-k+1 


25 ales 


Thus any number of ABM's equal to or less than (Xy) wax is acceptable. 
If the ratio X/T. satisfies the following inequality, the solution 


given in this paper can be used. 
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N-l i 
< i x 
My Ss tena ! 
i=l k=l 
For the example of section 7 
X/T, = 340/150 = 2.3 


and 


3 amet: 
P+ © tm (0.2 + 0.8p, )) = 1+ 0.92 + 0.74 + 0.53 
i=l k=l 


tt 1 it tt 


Dasadeye 
Since 2.3 < 3.2, the solution presented in this paper is applicable. 

9. Single Launchings Versus Multiple Launchings 

One of the assumptions made in this study was that only one 

ABM would be launched in defense of an undestroyed silo on any parti- 
cular stage, and that multiple launchings of ABM's would never be 
preferred over single launchings. The term "multiple launchings" 
used throughout this paper means the use of more than one ABM against 
a single RV aimed at an undestroyed silo on a given stage. The 
validity of this assumption is dependent upon the values of the para- 
meters of the problem. The purpose of this section is to devise a 
method for determining when the assumption is valid. Multiple 
launchings of three or more ABM's per silo per stage will not be 
considered, since this study is only concerned with those cases 
where single launchings are optimal. It follows that if single 
launchings are preferred over double launchings, then they are 
preferred over any other type of multiple launch. 

Consider eq. (3.40). An increase in x causes an increase in 


Tt given by 


Bo 


ap_(B + (1-q)""*) 


AX, (3.41) 


ee 
Ba a Aqp N 


A question one might ask is, given one additional ABM, should it be 
used according to the solution algorithm presented in this paper or 
should it be used to increase the defense of an undestroyed silo at 
stage one? To answer this question it is necessary to compare the 
marginal gain in Ts resulting from both options. 
If AX = 1, then ATS is equal to the coefficient of AXy in eq. 
(3.41) where AT* is the increase in the objective function if a new 
optimal solution (for single launchings) is computed. In the case 
where the additional ABM is used on stage one, the increase in Ts 
(call it D) is equal to the probability that a particular undestroyed 
silo survives stage one given that two ABM's are launched in its 


defense minus the probability that the same silo survives stage one ~ 


given only one ABM is launched in its defense. Thus 


(1-4 (-p,)”)- (2-4 Gp) 


D = qp,(1-p,). 


D 
(3.42) 


Consider the case where n=l, that is, where defensive firing begins 


on stage one. Using eq. (3.41) 


2 
qp, + 4Py qp, (1 a) 


i 
O 1+ qp, 1+ qp, 
Also D = qp, (1 - Py): 


Note that AT* > D 
since (1 + q)/(1 + qP,) > 1 and (1 - Pps Lea woieh 1S ‘o.say that” for 
Py > 0, single launchings are always preferred over double launchings 


on stage one given that defense begins on stage one. 
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Next consider the case where n=2. Then 


2 
d PoP; t+ aCl-q)p, dpy(1 - 4 + apy) 


AT* = = 
O 1+ 4P5 Lag 4P5 


D does not change since the same change in Xy is being considered 
throughout. Thus 
D = qp,(1-p,). 
The assumption of single launchings is valid in this case so long as 
po (1 - q + qp,) 
~I*q, = p, (1-p,)-. 

Therefore, for any specific set of parameters the following 
general comparison may be made to determine if the set of parameter 
values are such as to justify the use of the solution presented in 
this paper. If 

AT* > ap, (1-p,) 
for that value of n as determined by step one of the solution algorithm, 
then the initial assumption denying multiple launchings is valid on 
probabilistic grounds. There may be physical, engineering, or other 
restrictions which prohibit multiple launchings. 

Using the parameter values as given for the example of section 7, 
Where n=4 

ATS = 0.104 and Hb = w.192 
which indicates that at least one of the ABM's designated for stage 
four should be used on stage one if there is no valid reason to pro- 
hibit multiple launchings. Thus, the example, although sufficient 
for displaying the use of the solution algortihm, does not produce a 


maximum survivability if multiple launchings are allowed. Therefore, 


oy 


from a probabilistic point of, view, there are cases where multiple 
launchings are desirable even though some silos are undefended in 


the early stages of the attack. 
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IV. EXTENSIONS OF THE PROBLEM 


A. MULTIVALUED TARGETS 

Since this section involves equations which are analogous to model 
II, and because a detailed analysis of this subject is not intended, 
only those equations which characterize the problem and those comments 
which could be useful for obtaining a solution are given. The next 
section assumes silos of three different values for discussion purposes 
but any number of values (up to Ty) is possible. 

Consider each silo ot the farm to have associated with it a value, 
Va> Yaz OF Vary where, without loss of generality, it can be assumed 
that 


Such values can be related to the retalitory capabilities of the ICBM's 
located in the silos. 


There are now three distinct classes of silos to be considered, 


with es Tae and ue numbers of silos in each class. The total 


number of silos in the farm is just the sum of the numbers of silos 


in each class, 


cnt LT eel 
1 ee a le 


All other assumptions of the basic problem are valid in this generalized 
case. Let (TW) 5 be defined as the total worth of all undestroyed 
silos immediately prior to stage (N-i). Then 

I Il ie 


(TW) = Vary = Vern a Vrrtn 


a, 


1. The Objective Function 
The purpose ‘of the ABM system is to defend the farm so as to 
maximize the capability for retaliation. Thus, in a development 
analogous to the development of eq. (3.12), the generalized objective 
function is 
II 
y 


a N_m 2 eae 
max (TW) | = . ve | (1-4) T, + q > (l-q) pd (el ) 


j=l 


where a ds and ar are the decision variables representing the 


numbers of ABM's to be launched in defense of each of the three classes 
of silos on stage j. Such an objective function maximizes the total 
capability for retaliation following an attack. 
2. The Constraints 
The generalized constraints can be written directly from eq. 
(3.16) of model II. They are 


ITi N-1 


m 
Ni, aN 
m=I i=0 
m : j-(N-itl1)  m m i 
(2) d" -q E (1a) p,d, S Ty (1-a) (4.2) 
j=N-it1 
i=l, tae 
m=I, III, III 
(3) bee 0 ier Rie 
m=f, II, III 
3. Analysis 


Note that, in eq. (4.1), the quantity 


N IIl 
(1-q) ye Vat 


m= 


m 
N 
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is fixed for a given set of parameters. Thus, optimization (maximization 
of (TW) .) shall be concerned only with the remaining terms. Let qe be 

J 
the coefficient of qi in each of the remaining terms, i.e. 


Gir =v (l- ie is 
j a q) P, 


As in model II, an evaluation of the C. factors leads to the optimal 
solution provided 


> - JoeD) . 
ee (l-q) (see section III.B.5) 


One useful way to compare the cits is in the form of a matrix 
where of = c.. Each row of the matrix relates to a stage of the 
attack and each column is associated with a distinct class of silos. 
Thus, the matrix for this section is N x 3 whereas the matrix for 
model II is an N-dimensional vector. 

Assume that each element of the matrix is different. Then there 
exists an ordered sequence of matrix elements such that the first 
member of the sequence is the largest element of the matrix, the 
second member of the sequence is the second largest element of the 
matrix, and so on. With the earlier assumption that Vy Se Vat < Vaqr? 
the first ABM's would be assigned to defend silos of class LII on stage 
one. The analysis could continue as in model II, but a separate 
algorithm would be necessary for every feasible sequence of elements. 
The number of feasible sequences increases rapidly as the number of 
classes of silos increases and the method presented in this paper 
becomes inefficient. A more general approach is necessary to obtain 


an efficient solution to this generalized problem. Obtaining such a 


solution is recommended as a topic for further study. 
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B. OTHER TARGETS 

It was assumed throughout this paper that the ABM's and associated 
defensive equipment such as the radars and computers were immune to 
attack by enemy RV's. Of course, in a real situation this probably 
would not be true. The enemy might decide it would be to his advantage 
to strike the defensive system first and then use fewer RV's on the 
silos. In such a case a much larger problem exists than was presented 
earlier. It would not be sufficient to merely assign radars, computers, 
and ABM's values as eapeaEemsdinde the values of radars, computers, and 
ABM's are all interdependent. If no ABM's are available, radars and 
computers are useless, and if the radars or computers are destroyed, 
unlaunched ABM's are also of no value. Also, a simple, saturation, 
offensive strategy could not be assumed for in this case there exists 
a trade-off between attacking silos and attacking the ABM system in 
terms of motel silo destruction. It is only natural to assume that 
the enemy would analyze this trade-off and attack in such a manner as 
to maximize total silo destruction. Thus, for useful results, it 
would be necessary to peeinize offensive and defensive strategies, 
maleate oust 

There are numerous other variations to the basic problem investi- 
gated by this paper. There are also numerous other methodologies 
available ane theory, decision theory, computer simulation) for 
formulating, Baten and solving the problem variations. The ABM 
question is a Foret le field for analysis and is sure to receive a 
great amount of attention in the future. This paper has been a modest 
attempt at solving a simple form of a complex problem, and has tried to 


present some insight into the nature of the greater problem. 
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APPENDIX A 


THE ANALYTICAL DIFFICULTY OF FRACTIONAL ABM'S 
Let P represent the probability that a target survives three 
consecutive, independent, attacking RV's given that a single ABM 
is launched in defense. Then 


P = (1-q) (1-q) E = (1-p)| 


where q and p are defined as in part B, chapter II. 
let R represent the probability that a target survives three 
consecutive, independent, attacking RV's given that 1/3 of an ABM 


is launched at each RV. Then 
ee E -@ (1-p/3) °. 


Now the problem remains to determine whether P > R or R> P. Assume 


R> P. Then 
1 Z 
=" quel 39) > (l-q) | 1 - q (Lop). 
Expanding both sides and collecting terms yields 
2 S| 2 Cae Zed 
LHS = 1 - 3q + 3q° - q+ pa + pq - 2pq’. +p q/3 - pq/3 
3 
40: oe 


RHS 


D 3 8 2 
ieee dete5g - q &£ pq pq - 2pq 


The first seven terms of the LHS cancel the seven terms of the RHS, 


leaving 
sachs - SOE + p aw27 > OF 
But pq-/3 > p-q?/3 
since w= gq < lt 
and per /27 0. 
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Thus R = P only if p = 0 and q = 1. Such values taken together are 
totally outside the realm of interest. Thus it can be concluded, 
for p = Oymand adinyeoe= cn is 
| Roe. 
So although it is mathematically convenient to do calculations 
using fractional ABM's, it may well lead to inflated values for 


target survivability. 
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APPENDIX B 
DETAILS IN THE DEVELOPMENT OF MODEL I 


A. REDEFINING THE OBJECTIVE FUNCTION 


From eq. (3.1) 


N 
Z = YY: (B.1) 


where Yow = gi...- qped. (eq. (2.4)). 
J J J 3 


Therefore Y, = qT, - qp,d, 


Yo = 4To - qPod, 


yoy iy Oe en) 

N 
and z= y rT. = p.d.). B.3 
oe p,d,) (B.3) 

ml 

However, from Lemma 1 

N-j 2 etal 
TaeT (1-g) 1 +4 5 (-g) ‘1p de. (B.4) 
J k=jt+1 


So (B.3) becomes 


N : 
zZz=q>b ne (1- op ad J + q : (1- aay Sen) 4, | -q . p.d 


j=l k=j+1 jel J J 
N N N N 
2 -(jtl 
Zz = qty x (l- ae ee ag y (1-q)* Gj ody q 2 p.d. 
ie jal k=jt+l j=l JJ 


N 
N- k-(4+1 
Ze que. = (l=q) el q : | ; (I=q) (3 dy : pd, | (B. 5) 


j=l j=l ~ k=jtl 
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Letting N-j=i, the first term can be written 
N-1 


pee Clas 
i=0 


N 
which is identical to the expression for dD, ineqewse5).. Lt now 

remains to separate the decision variables. Note that d, is wii Gon ly 
one term of the expansion, d, is in two terms, d, is in three terms, 


and so on. 


Expanding the bracketed part of eq. (B.5) yields, for 


en mL 2 __\N-2 

j=l: dpod, + aQ1-q) p3dg + qQl-q) pyd, +... + q(l-q) pydy - pydy 

j=2: (1- ve d, + q(1- e d, + + q(1- ye d_ - pad 

jee: q(I-q) p,d, + q(l-q) p,d, + ... + q(1-q Pyty 7 Pods 

as _\9 a Nae : 

io: q(1-q) py dy, + ... + q(l-q) Prytig pd. 

j=N-1: q(1-q)"p a =p ‘ 
NN N-1 N-1 

j=N: d 


~ Py“y 


Summing terms common to each oe and multiplying through by q produces 


; s 


l N-2 
t t 
-9p,;d, - dp (l-a)d, + dp, g ~ (l-q) - 1 | da +... + apy E zy (l-} -1 a. 
C= E=C nN 


0 
(B.6) 


The first term above;is;D, as used in eq. (3.5). The general expression 


if 


d_. cam Be written 


for any of the terms above involving dys sees N 


f% k-2 . 

ap, (4 C-a)®- 1] 4, 

t=0 
which is equivalent to the general term in eq. (3.5) when the sub- 
stitutions t=j and k=N-i are used. The proof of eq. (B.6) will not 
be given. It is analogous to the.inductive proofs used for Lemmas l 


and 32 
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In conclusion, the objective function can be written 


N 
z= %)} Y. or equivalently, 
ie 


N 
i 


D, st: Did, Sg Dod, heme 4 et Dy dy 


with the coefficients of the Slegmit as defined above. 
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